Fast Zonal Field Dynamo in Collisionless Kinetic Alfven Wave Turbulence 
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The possibility of fast dynamo action by collisionless kinetic Alfven Wave turbulence is demon- 
strated. The irreversibility necessary to lock in the generated field is provided by electron Landau 
damping, so the induced electric field does not vanish with resistivity. Mechanisms for self-regulation 
of the system and the relation of these results to the theory of alpha quenching are discussed. The 
dynamo-generated fields have symmetry like to that of zonal flows, and thus are termed zonal fields. 
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The dynamo problem, i.e. the problem of understanding the origin and generation mechanism for astrophysical and 
geophysical magnetic fields, remains one of the major unsolved problems in classical physics Until recently, the 
kinematic quasilinear theory of the alpha effect provided an attractive framework for formulating a solution to this 
!_} 1 problem In this approach, quasilinear iteration was used to relate the turbulence-induced mean EMF, i.e. (?x B) 
to a(B_) — f3( J), where a, the dynamo drive term is in turn proportional to the helicity of the fluid turbulence, i.e. 
(v-uj). Recent computational and theoretical research has, however, revealed that the kinematic theory is invalid, since 
the small scale magnetic field grows on fast, inertial range time scales, and so cannot be neglected when considering 
the evolution of the mean-field ||. The 'bottom line' of this research is that an , the kinematic alpha coefficient 
discussed above, is substantially reduced by self-consistent small scale field effects, so that the net alpha effect is [|J 
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Here vao = (M}/V^Po is the Alfven speed in the mean field, (v 2 ) is the mean square turbulence velocity, R m is 
the magnetic Reynolds number and ax = ^2(u 1 i±L)kT c ,k- T c ,k is an eddy self-correlation time. Eqn. (|l]) states 
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that the alpha effect is quenched for va > (1/ WRm)vrms, i-e. for very small (B), since R m >> I in most relevant 
applications. Note that the above expression for a may be re- written as a — aK'V/iv + T c v Ao\i where r c is the 
integral scale turbulence correlation time. This expression suggests that one way of interpreting the alpha 'quench' 
•£3 ■ is that despite naive expectations of turbulent mixing, etc., the collisional resistivity ultimately controls the rate at 
which 'stretched and twisted' fields are folded to increase (£?). This constraint necessarily precludes the possibility 
• • . of a fast dynamo, i.e. one which operates on time scales independent of r\. Another possible interpretation of the 
, £^ ' quench is that of 'Alfvenization', whereby the growing (B) converts eddys into Alfven waves (which are intrinsically 
non-kinematic), thus initiating the quench. Indeed, with this in mind, it is interesting to note that the EMF for an 
individual, circularly polarized visco-resistive MHD Alfven wave is < v_ x B_ >= (i]/2B )(l — v/r/)k 2 k z A 2 ||. Apart 
from the 'cross-helicity' factor (1 — v jr\) and a somewhat different expression for the helicity spectrum, this formula 
is identical to that given by Eqn. (Q), thus re-inforcing the notion that Alfvenization is at work in the alpha quench. 

In this paper, we report on the theory of a fast, collisionless dynamo. In view of the above discussion of Alfvenization, 
we proceed to directly consider a dynamo driven by Kinetic Shear Alfven Wave (KSAW) turbulence ||. Here, 
the dynamo instability mechanism is that of a modulational instability of the KSAW spectrum. In this case, the 
dynamically generated field is orthogonal to an externally prescribed mean field B_ a = BqZ, along which the KSAW's 
propagate Q| . Collisionless dissipation, via Landau damping, locks in the dynamo-generated field at a rate independent 
of 77, thus facilitating fast dynamo action. Since, assuming radially inhomogeneous turbulence, the dynamo-generated 
fields have azimutal symmetry and are thus analogous to zonal flows ||, we hereafter refer to them as zonal fields. 
As the current is carried by the electrons in KSAW turbulence, we start from the drift-kinetic equation 
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Here E\\ = — V||0 — (l/c)dA/dt, where <j> is the electrostatic potential, A is the z component of the vector potential, 
and V|| = d/dz + ¥.(A/B ) x z ■ V. Neglecting electron inertia, averaging Eqn.(j|), and using Ampere's Law then 
yields the mean field Ohm's Law for collisionless plasma: 
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Terms (E z n/n ) and dT/dr refer to parallel electron acceleration and the spatial transport of parallel current, respec- 
tively. The flux r is given by (l/fi e ) J dvuvu ((Ve{4> — (v\\/c)A\\))f). Proceeding as in quasilinear theory, we can use 
the coherent response 
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to calculate (E z n/n ) and T, so that 
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Here f/ ; = uA\\/ck z and <^> and t/» are related by ipk — (1 + kj_p^)(j)k- Only the non-adiabatic piece of / contributes to 
(E z n/n ) and T. Note that the flux of parallel current is proportional to a spectrally averaged factor of (k ± x z)/fc z , 
which requires that the turbulence have a net spectral chirality, in order for T to be non-zero. This property is the 
manifestation of helicity in the KSAW dynamo problem. Note also that T is independent of resistivity ry, since KSAW 
Landau damping now provides the irreversibility which permits 'fast' transport of current! 

The stability of the KSAW spectrum to a zonal magnetic field perturbation SBg = 5Bo(r)9 may be determined by 
modulating Eqn. i.e. 



1 d r; .. 5(E z nln ) . 5T r 



5(A) 



'5(A) 



(5a) 



where 5(A) is the associated modulation in the vector potential, i.e. 5Bg(r) = —85(A) /dr. Anticipating the use of 
methods from adiabatic theory, (E z n/n ) and T r are conveniently re-expressed in terms of the KSAW action density 
Nk as 
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where — Wk/^k and the KSAW energy density (normalized to the thermal energy density n a T e ) is 
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Using the relation between potentials 4>k and Au/., this can be re-expressed as: 
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Here p s = c s /Qi and c 2 = T e /rrii. In this formulation, then, the modulations 8(E z n/n ) and <5r can now be calculated 
simply by computing the modulational response SNk/S(A). 

The modulation in the number of KSAWs induced by the dynamo generated zonal field 8(A), i.e. 5N k /S(A), may 
be obtained by using the wave kinetic equation to determine the response SN k to ^(^4). The action density of a wave 
packet evolves according to the wave kinetic equation 

— + v g ■ VN k - Vuj k ■ — = ~j k N k , (6) 

where jk is the wave damping decrement. The dynamo generated field refracts the underlying KSAWs and thus induces 
spectral modulations. In the presence of dynamo-generated zonal fields Uk — > ^ k + 6ui k . Here uj k = k 2 v\{\ + k\p 2 ) 
is the linear wave frequency and 

Sujk = r~| k g SBe, (7) 

fez-Do 

is the perturbation in the wave frequency induced by the zonal dynamo field. Linearizing Eqn. (^) then gives the 
action density modulation 
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Here q is the radial wave number of the modulation 6(A) and fL is its frequency. (Nk) is the mean KSAW action 
distribution. Eqn. (||) may then be substituted into Eqns. ( papqfk| ) to obtain the general expression for fl q , which 
is (neglecting resistive dissipation) 
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Here d e = c/oj p . e , p e is the electron gyroradius and f = f a / v the- Now, noting that the radial group velocity is 
Vg.r — duJk/dk r — (k 2 v\/oJk)k r p 2 and considering the limit qv g ^ r > 7, il q (note here f2 g is small, i.e. of 0(E 2 ), and 
qv gr > 7 is consistent with the assumption of weak wave damping) , we ultimately obtain the zonal field eigenfrequency 
flq, i.e. 
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Several aspects of this result are of interest. First, note that zonal field growth requires a normal population profile, 
i.e. d(Nk) / dk 2 < 0, similarly to the case of zonal flow generation by drift waves. This condition is virtually always 
satisfied for Alfvenic MHD turbulence. Second, note that zonal field growth ImQ, q is independent of 77, since electron 
Landau damping provides the requisite dissipation. However, concommitant with this is the fact that the dynamo 
drive is ultimately proportional to \E z k\ 2 /k 2 — \{<f> — ip)k\ 2 , and thus depends directly upon field-fluid decoupling 
due to finite gyro-radius (and possibly finite ion inertial layer width, i.e. finite k 2 c 2 /oj 2 ^ effects. This is a direct 
consequence of the fact that finite En is required for Landau damping. Since, in turn, Landau damping provides the 
necessary mechanism for irreversibility, the KSAW dynamo necessarily then requires non-zero Eh on KSAW scales. 
Note also that while a mean T (and thus a mean field KSAW dynamo) requires a net chirality or helicity in the KSAW 
turbulence, zonal field growth does not. This is because zonal field formation is a spontaneous symmetry breaking 
phenomenon, whereby the system acts to reinforce an initial 8B§ of either sign, but does not amplify the total flux. 
In this respect, the zonal field amplification process is more like the small-scale dynamo than the mean field dynamo. 
However, we hasten to add that a broad spectrum of growing zonal fields, with (pi/L±) < qpi < 1, can be expected 
(N.B. Here L± is the scale of the Alfven wave spectrum inhomogeneity in the direction perpendicular to B z). Thus, 
the zonal field dynamo field is not restricted to hyper- fine scales such as qpi ~ 1, etc. A further observation is that 
the zonal field dynamo mode are predicted to oscillate as well as grow, i.e. Q q is complex. A net chirality (i.e. finite 
spectrum averaged kgk z ) is required for Re fl q y^0. 



3 



Since KSAWs (with finite k r ) will be refracted by zonal fields, the KSAW spectrum will necessarily be modified 
by the growth of the dynamo-generated fields. Thus, a self-consisted KSAW dynamo theory, which treats both the 
zonal field generation as well as their back-reaction on the KSAWs, is called for. The effects of random refraction of 
KSAWs by the generated fields can be described by the quasilinear wave kinetic equation 

d{N k ) d duj k *t \ /at \ /m •> 

~ar ~ = (10a) 

which, using uj k and N k as given by Eqns. (@,|[), can be simplified to 

d(N k ) d d(N k ) 
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R(Q q - qvg r ) = Jk/[(Cl g - qvg r ) 2 + t|] = 7fc/ (qv gr ) 2 . (lOd) 

Here, the irreversibility intrinsic to D k is provided by KSAW electron Landau damping via 7&. Since zonal field 
frequencies are low (i.e. Vl q ~ (SB / B ) 2 ), D k is effectively non-resonant in character. Note the effect of random 
refraction is to 'diffuse' k r to higher values. This refraction occurs as KSAW packets traverse the layered structures 
of zonal fields. Since wave dissipation is likely to be stronger at high k r , this constitutes a route for feedback of the 
fields on the KSAW intensity, as well as its spectral distribution. 

Having obtained the evolution equation for (Nk), we can now self-consistently describe the coupled evolution of 
the dynamo driven zonal field spectrum and the KSAWs. These may be described by coupled predator-prey type 
equations for \SBg q \ 2 and (Nk), which are: 

^\SB 0>q \ 2 =2q 2 (j g -r 1 ), (11a) 



and 
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where 
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Here 7 ? = 7 g ((A)) is derived directly from Eqn. ( j9b| ) for Im(Cl q ). Note that Eqn. (11a) implies that a critical KSAW 
intensity level (N) cr u ~ r) is required for zonal flow growth via modulational instability. For weak 77, this implies that 
only a very modest level of KSAW excitation is required for the dynamo. Moreover, in regimes of weak dissipation, 
it suggests that the 'marginal point' for the KSAW dynamo will scale with 77, despite the fact that 77 is otherwise 



completely irrelevant to the dynamics of zonal field generation. Eqn. (lit) suggests that as zonal fields grow, high 
k r components in the KSAW spectrum will be generated, thus ultimately quenching KSAW energy via coupling to 
dissipation. Finally, it should be noted that the expression for j q given in Eqn. flllcj ) assumes qpt < 1. For qpi ~ 1, 
additional FLR factors enter which force 7 Q to decay for qpi > 1. 



A moment's consideration of Eqns. (11a, lib 11c) naturally begs the question of what happens in the limit of 
X] — * since, in this case, there is no apparent control on zonal flow growth. We speculate here that as r\ — * 0, zonal 
fields may be subject to collisionless reconnection instabilities (i.e. tearing modes), which limit zonal field growth 
without requiring resistive dissipation. Such instabilities are analogous to Kelvin-Helmholtz type instabilities which 
may limit zonal flow growth in the zero collisionality limit || . Alternatively, the zonal fields may regulate themselves 
via feedback on d(N k ) / dk 2 (i.e. by modifying the wave spectrum |l(J) or by trapping of KSAW packets |11|. The 
detailed dynamics of these collisionless reconnection instabilities and of the system's behavior for 7/ — > will be 
discussed in detail in a future publication. 



4 



In conclusion, we have demonstrated that a zonal field dynamo can be driven by collisionless kinetic Alfven wave 
turbulence. This dynamo is 'fast' in that the rate of requisite magnetic reconnection is determined by electron 
Landau damping and thus is independent of the collisional resistivity. Concommitant with this, non-zero E\\ (i.e. as 
due to FLR or finite ion inertial layer width) is required for magnetic field amplification. The zonal field dynamo 
self- regulates via refraction-induced diffusion of the KSAW spectrum toward high k r . These predictions should be 
amenable to investigation in laboratory experiments. Finally, it is interesting to note that for typical interstellar 
medium parameters (T e ~lev, B o ~l0~ e G, L ~ 1 pa, etc.) and taking B / B o ~10~ 3 , Eqn. ( pb| ) predicts that the 
magnetic field grows at a rate 7 o ~10~ 17 sec~ 1 , consistent with that needed to achieve equipartition in 10 9 years. 
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